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First MINOS Results from the NuMI Beam 

Nathaniel Tagg, for the MINOS Collaboration 

Tufts University, 4 Colby Street, Medford, MA, USA 02155 

As of December 2005, the MINOS long-baseline neutrino oscillation experiment collected data 
with an exposure of 0.93 x 10 20 protons on target. Preliminary analysis of these data reveals a 
result inconsistent with a no-oscillation hypothesis at level of 5.8 sigma. The data are consistent 
with neutrino oscillations reported by Super-Kamiokande and K2K, with best fit parameters of 
Am| 3 = 3.05±g;|§ x HT3 and sin 2 2(9 23 = 0.88±^ 2 . 
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I. INTRODUCTION 

The MINOS long-baseline neutrino oscillation ex- 
periment was designed to accurately measure neu- 
trino oscillation parameters by looking for Vn disap- 
pearance. MINOS will improve the measurements 
of Am 2 23 first performed by the Super-Kamiokande 
dill and K2K experiments 0. In addition, MINOS 
is capable of searching for sub-dominant — > v e os- 
cillations, can look for CPT-violating modes by com- 
paring to oscillations, and is used to observe 
atmospheric neutrinos [5j. 

The MINOS experiment uses a beam of created 
at Fermilab National Laboratory and directed at the 
Soudan mine in northern Minnesota, a distance of 
735 km. The composition and energy spectrum of the 
beam is measured in two detectors, the Near (1 km 
downstream from the target) and the Far (735 km 
downstream), allowing for precision measurements of 
the spectral distortion of the beam. 



II. THE NUMI NEUTRINO BEAM 

The Fermilab Main Injector has a minimum cycle 
time of 1.87 s, with a maximum intensity of 4 x 10 13 
protons per pulse for a maximum average power of 
0.4 MW of power on the target. Protons are extracted 
from the Main Injector in a single turn, taking ~ 10/is. 

The NuMI neutrino beam |6J is created by directing 
these protons onto a water-cooled segmented graphite 
target. Secondary pions from the proton-carbon in- 
teractions are deflected into the forward direction by 
two parabolic focusing elements (horns) before being 
allowed to decay in a 675 m evacuated decay volume. 
Undecayed secondaries are stopped by an absorber 
wall at the end of this volume. The primary proton 
beam is monitored for position and intensity on target. 
The secondary hadrons and muons are monitored for 
position and intensity to ensure good alignment and 
composition of the beam. 

The NuMI beam has been in operation since late 
2004. Near the end of 2005, a total exposure of 
0.93 x 10 20 protons were delivered to the target. This 
exposure was used in the following analysis. At the 
end of the run period in March 2006, the maximum in- 



tensity delivered to the target was in excess of 25 x 10 12 
protons per pulse, with a maximum target power of 
250 kW. 



III. THE MINOS DETECTORS 

The MINOS Near and Far detectors are constructed 
to have nearly identical composition and cross-section. 
The detectors consist of sandwiches of 2.54 cm thick 
steel and 1 cm thick plastic scintillator, hung verti- 
cally. The polyethylene scintillator is in the form of 
4 cm-wide strips, with a co-extruded Ti02 reflective 
coating. Along one side of each strip a groove holds 
a glued wavelength-shifting optical fibre. Scintillation 
light created in the scintillator is caught by the fibre, 
shifted, and transported efficiently to the end of the 
strip. Clear readout fibres carry the light to multi- 
anode photomultipliers for readout. 

The detectors act as tracking, sampling calorime- 
ters. Strips in adjacent planes are oriented orthogo- 
nally, allowing events to be reconstructed in two trans- 
verse views. Both detectors are equipped with magnet 
coils which generate ~1.2 T toroidal magnetic fields, 
which act to contain long muon tracks and provide 
curvature information for estimating energies. 

The Near detector has a total mass of one kiloton, 
and 282 of these steel planes, 153 of which are instru- 
mented. It uses sampling electronics to distinguish 
neutrino events in time, due to the large instantaneous 
intensity during a beam spill. The Far detector masses 
5.4 kilotons, and consists of 485 steel planes, 484 of 
which are instrumented. Because of the low rate in 
the Far detector, the strips are read out via an 8- fold 
optical multiplexing. 



IV. DATA SELECTION 

Charged-current events in the MINOS detectors 
appear as long muon tracks accompanied by short 
hadronic showers near the event vertex. The two 
other distinguishable event classes are neutral-current 
events, which appear only as short, sparse hadronic 
showers, and v e charged-current events, which appear 
as short, dense electromagnetic showers. For this anal- 
ysis, only charged-current events are considered. 
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Neutrino events are selected by first taking time- 
coincidence with the beam spill. This is performed 
by hardware trigger in the Near detector. In the 
Far detector candidate events are required to occur 
at 2449 ± 50/is after beam spill (the time-of-flight 
for neutrinos traveling to Soudan). Charged-current 
i/fj, events require a well-reconstructed track in both 
scintillator views, and confinement of the track ver- 
tex to the defined fiducial volumes of each detector. 
The curvature of the track is selected to be consis- 
tent with a ji~ . Finally, charged-current events 
must pass a particle identification cut, which relies 
on event length, the proportion of calorimetric energy 
in the shower, and the track dE/dx. These variables 
are assembled into a probability density functions and 
used to create a likelihood-based particle ID param- 
eter shown in Figure ^ A cut is performed on this 
variable to select a pure CC sample. The total ef- 
ficiency and purity of this selection is shown in Figure 
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FIG. 1: The Particle ID Cut. The solid line shows the dis- 
tribution of Monte Carlo charged-current events and the 
dashed line shows the distribution of Monte Carlo neutral- 
current events as a function of the particle ID parameter 
used for the selection. The cut value is shown. 

For selected events, the neutrino energy is recon- 
structed by taking the sum of the muon energy and the 
shower energy. Muon energy is found using the range 
of the muon (if contained) with an accuracy of 6%, 
or by using the muon curvature (if uncontained) with 
an accuracy of 10%. The shower energy is found us- 
ing calorimetry, with an uncertainty of approximately 
55%/VE. 



V. THE NEAR DETECTOR SPECTRUM 

The Near detector neutrino sample is used as the 
control in the experiment, providing the unoscillated 
beam spectrum. One of the advantages of the NuMI 
beam is that the target may be moved relative to the 
horns, allowing a varying focus of the pions and there- 
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FIG. 2: Purity and Efficiency of the CC f M Event Se- 
lection. Efficiency and purity are shown as a function of 
reconstructed event energy, using Monte Carlo. Contami- 
nation events are from NC events, almost entirely at low 
reconstructed energy. 



fore different neutrino energies. Three such settings 
are used to create 'high', 'medium' and 'low' energy 
beams. Although the bulk of data was taken in the low 
energy configuration (to maximize sensitivity at the 
oscillation minimum) data was also taken in the other 
configurations to test beam modeling in the Near de- 
tector. 

Monte Carlo simulation of the beam predicted the 
neutrino flux for all beam configurations to within 
10%. The residual disagreement was a strong func- 
tion of beam configuration rather than neutrino en- 
ergy, implying that the simulation error was due to 
the modeling of the secondary pions. An empirical 
function was used to modify the secondary pion simu- 
lation as a function of Xf and p t . This function was fit 
to data from all configurations. These data and their 
associated Monte Carlo spectra are shown in Figure 

El 

The measured Near detector spectrum was used to 
predict the Far detector spectrum. First, the observed 
Near detector reconstructed energy spectrum was cor- 
rected for purity and efficiency to obtain the true neu- 
trino spectrum. 

Second, the spectrum was corrected for the differ- 
ence in solid angles subtended by the Near and Far 
detectors. The smaller range of angles observed by 
the Far detector results in a narrower spectrum rela- 
tive to the Near detector. This difference can be ro- 
bustly estimated by using knowledge of the pion two- 
body decay kinematics and the geometry of the NuMI 
beamline. Using the beam Monte Carlo, a matrix is 
created that transforms a true Near detector spectrum 
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FIG. 3: Near Detector Neutrino Energy Spectra. The 
spectrum of observed charged current dots. Hatched areas 
indicate the Monte Carlo prediction of the spectrum, after 
tuning. The predicted neutral current contamination is 
shown in blue. The top plot shows the data from the 
low energy beam configuration, the middle plot shows the 
medium energy configuration, and the lower plot shows 
the high energy configuration. 



VI. FAR DETECTOR RESULTS 

The observed reconstructed neutrino spectrum at 
the Far detector is shown in Figure 0] The data do 
not agree with the unoscillated prediction; this null 
hypothesis is excluded at 5.8 sigma. When fitted to 
a two-neutrino oscillation hypothesis, a fit is found 
with reasonable probability within a parameter space 
shown in Figure|SJ In this fitting procedure, only sta- 
tistical errors were considered. 

Systematic errors were estimated by taking Monte 
Carlo data samples simulated near the best fit point 
for the data. Systematic errors were introduced to the 
Monte Carlo and not corrected in analysis. The result- 
ing shift in the best fit point is taken as the systematic 
error on that parameter. The systematic errors con- 
sidered are shown in Table [I] These systematics are 
briefly discussed below. 

Normalization changes the overall expected num- 
ber of events at the Far detector, based upon a 
combination of uncertainty on protons-on-target 
and the fiducial mass. This is known to 4%. 

Muon energy scale is the uncertainty on the aver- 
age energy of muons due to incorrect steel den- 
sity or magnetic field modeling. This is known 
to 2%. 

Relative shower energy scale is the degree of con- 
fidence to which the hadronic calorimetry in the 
Near detector matches that in Far. The energy 
scale between these two detectors after calibra- 
tion should agree to within 3%. 

NC contamination is the uncertainty on the num- 
ber of neutral current background events erro- 
neously accepted into the charged current sam- 
ple. A conservative error of 30% is used. 

CC Cross section is allowed to vary in several ways: 
the KNO parameters governing the relative frac- 
tion of resonance to deep inelastic scattering 
fraction were allowed to change by 20%, and the 
axial mass for quasi-elastic events and resonance 
events were allowed to change by 5%. 

Beam uncertainty is the difference in measured pa- 
rameters given that the data is fit without re- 
tuning the Monte Carlo beam simulation to the 
Near detector data. That is, an incorrect beam 
matrix was used to extrapolate the beam to find 
the shift in fit values. 



into a true Far detector spectrum. 

Last, the predicted true Far detector spectrum is 
simulated by Monte Carlo to create an expected recon- 
structed neutrino spectrum which can be compared to 
data. 



Intra-nuclear rescattering gives the uncertainty 
due to poor modeling of the energy loss of pi- 
ons scattered or re-absorbed as they leave the 
target nucleus. This energy cannot be detected, 
and so an uncertainty on this interaction creates 
an uncertainty on the absolute hadronic shower 
energy scale at both detectors. 
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result will significantly improve the world average on 

X 2 / n.d.f = 20.5/13.0= 1.6 



TABLE I: MINOS PRELIMINARY Systematic Uncer- 
tainties. Systematic effect on the data are shown together 
with the shift incurred by not correcting for these system- 
atics on the best-fit point of Am| 3 and sin 2 2#23 • 
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Table |U also shows the sum of these systematic er- 
rors taken in quadrature, as compared to the statisti- 
cal errors. It can be seen that for this limited data set, 
MINOS is statistically limited. Systematic errors have 
a negligible effect on the contour at this sensitivity. 
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FIG. 4: MINOS PRELIMINARY Far Detector Energy 
Spectrum. The predicted spectrum in absence of oscil- 
lations is shown by the solid curve. The dashed curve 
represents a fit to the data by an oscillation hypoth- 
esis with a best fit of Am! 3 = 3.05±o.55 x 10" 3 and 
sin 2 2(9 23 = 0.88±°;i 2 



VII. CONCLUSIONS 

MINOS has completed a preliminary analysis of the 
first ~ 10 20 protons on target, and found a result in- 
compatible with no oscillations at a significance level 
of 5.8 sigma. The signal is consistent with oscillations 
seen by Super-Kamiokande and K2K at the 90% con- 



FIG. 5: MINOS PRELIMINARY Allowed Region. The 
black solid (dashed) line shows the MINOS 90% (68%) 
confidence limit in Amis an d sin 2 2023- The 90% confi- 
dence regions from the Super-K zenith-angle analysis @], 
the Super-K L/E analysis Q and the K2K result Q are 
shown in red, green, and blue respectively. The contour is 
described by statistical errors only. 



Amp. 

MINOS currently has a data sample of approxi- 
mately 1.5 x 10 20 protons on target, and will be analyz- 
ing this expanded data set for summer 2006. MINOS 
intends to continue collecting data until at least 2010, 
at which time it's physics sensitivity is expected to be 
considerably improved. 
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